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ABSTRACT 

In this paper we show that a one-to-one correspondence exists between any dark 
energy model and an equivalent (from a cosmological point of view, in the absence of 
perturbations) quartessence model in which dark matter and dark energy are described 
by a single perfect fluid. We further show that if the density fluctuations are small, 
the evolution of the sound speed squared, c^, is fully coupled to the evolution of 
the scale factor and that the transition from the dark matter to the dark energy 
dominated epoch is faster (slower) than in a standard ACDM model if > (c^ < 0) . 
In particular, we show that the mapping of the simplest quintessence scenario with 
constant wq = Pq/pq into a unified dark energy model requires < at the present 
time (if wq > — 1) contrasting to the Chaplygin gas scenario where one has > 0. 
However, we show that non-linear effects severely complicate the analysis, in particular 
rendering linear results invalid even on large cosmological scales. Although a detailed 
analysis of non- linear effects requires solving the full Einstein field equations, some 
general properties can be understood in simple terms. In particular, we find that in 
the context of Chaplygin gas models the transition from the dark matter to the dark 
energy dominated era may be anticipated with respect to linear expectations leading 
to a background evolution similar to that of standard ACDM models. On the other 
hand, in models with > the expected transition from the decelerating to the 
accelerating phase may never happen. 
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1 INTRODUCTION 

A considerable effort has been devoted, in the past few years, 
to shape Unified Dark Energy (aka Quartessence) into a 
coherent phenomenological model. At the moment we are 
still lacking a solid theoretical motivation in favour of this 
unified scenario. Yet, the potential outcome of such models 
is enough to keep a significant amount of papers cyclicly 
emerging. The drive behind UDE is simple enough to be put 
into a few words: can dark matter and dark energy share a 
common unifying nature? Could they be different aspects of 
the same thing? Either way, an answer to these questions 
would tell us something of a fundamental nature. 

Evidence for the existence of a dark energy component 
has been steadily piling over the years. It is clear from ob- 
servations that most of the matter in the Universe is in a 
dark non-baryonic form (see, for instance, Tonry et al. 2003; 
Riess et al. 2004; Tegmark 2004; Spergel et al. 2006). How- 
ever, there is at present no direct detection of non-baryonic 
dark matter or dark energy, their existence merely inferred 
from their cosmological implications through gravitational 
effects. Hence, one should not discard the possibility of a 



single component simultaneously accounting for both dark 
matter and dark energy. 

Historically, the idea of UDE has sprung from the cos- 
mological properties of the Chaplygin gas (Kamenshchik, 
Moschella & Pasquier 2001; Bilic, Tapper & VioUier 2002; 
Bento, Bertolami & Sen 2002) , an exotic fluid with an equa- 
tion of state p — ~Ap^" (where A and a are positive con- 
stants). This special fluid has a dual behaviour: it mimics 
matter (p ~ 0) early in the history of the Universe and a cos- 
mological constant much later (with a smooth transition in 
between) which is highly suggestive of a unifled description 
of dark matter and dark energy. 

Although we could assume that the Chaplygin gas is 
nearly homogeneous in both the radiation and matter eras 
much like a standard variable lo-quintessence component 
(see for example Multamaki, Manera & Gaztanaga 2003), 
this would not explore the full potential of the Chaplygin 
gas as a UDE candidate. Hence, in the absence of any real 
clues regarding the microscopic properties of quartessence, 
we will make the simplest possible assumption: that it can 
be approximated by a perfect fluid whose properties are fully 
specifled once the equation of state is known. However, we 
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shall refrain from making any further assumptions regard- 
ing the form of its equation of state (other than those which 
are required by observations that is) and leave the discus- 
sion as general as possible (see, however, Bento, Bertolami 
& Sen 2004 for a different approach). 

An appealing feature of perfect fluid UDE models is the 
existence of formally equivalent models with a single scalar 
field, (j>, described by the action 



S = / d\^gp{X) 



(1) 



where p = p{X), p = 2Xdp/dX — p and — Vfi</>/\/2X 
are respectively the pressure, the energy density and the 
four velocity of the perfect fluid and X = V ^(l)^'^<f)/2. If 
p = x^^'^^^^^^'^^ the action describes a perfect fluid with 
equation of state p = wp with constant w ^ 0,-1 (for w = 1 
one obtains the usual massless scalar field action) while if 



p(X) = 



(2) 



one gets a generalized Born-Infeld action describing the dy- 
namics of the generalized Chaplygin gas. Here we are taking 
V i_i(j> to be timelike. In all these models a single real scalar 
field accounts for both dark matter and dark energy. 

The predictions of UDE models based on the Chaply- 
gin gas have been tested using observational data including 
high-z supernovae (Avelino et al. 2003a; Bean & Dore 2003; 
Bega et al. 2003; Colistete & Fabris 2005; Bento et al. 2004; 
Zhu 2004; Bertolami et al. 2005), lensing (Makler, Quinet 
de Oliveira & Waga 2003; Silva & Bertolami 2003; Dev, 
Jain & Alcaniz 2004), high precision CMB (Carturan & 
Finelh 2003; Bean & Dore 2003; Amendola et al. 2003; 
Bento, Bertolami & Sen 2003) and Large Scale Structure 
(Fabris, Goncalves & De Souza 2002a; Fabris, Goncalves & 
De Souza 2002b; Bean & Dore 2003; Bega et al. 2003; Sand- 
vik et al. 2004) . As we'll later see in more detail, the Chaply- 
gin gas attains very high sound speeds at recent times which 
has a significant negative impact on small scale structure for- 
mation. This was first studied by Sandvik et al. (2004) which 
concluded that in order to obtain the mass power spectra we 
observe today, the parameter a in the Chaplygin gas model 
had to be extremely fine tuned around zero (q = being 
the ACDM limit). However, it has been shown by Bega et 
al. (2003) that this problem could be alleviated (but not 
solved) by adding baryons to the mixture. 

Yet, there is a caveat to most of these results: linear 
theory has been taken as a good approximation when com- 
paring with cosmological observations. Ordinarily, this is a 
valid assumption to make; for example in standard ACDM 
models the non-linear collapse is not expected to have a 
large impact on the evolution of the average Universe. How- 
ever, as discussed previously by Avelino et al. (2004), the 
small scale structure of the Chaplygin gas may infiuence 
the background evolution of the Universe, rendering linear 
theory invalid even on large cosmological scales. Hence, in 
order to confront perfect fiuid UDE models with observa- 
tions, linear theory may not enough; we may have to solve 
the full Einstein field equations which is obviously an enor- 
mous task. Nevertheless, in this paper we will show that 
the overall effect of the non-linear terms on the background 
evolution of the Universe may be understood using simple 
considerations. 



2 LOCAL VS. GLOBAL EQUATIONS OF 
STATE 

The dynamics of a homogeneous and isotropic Friedmann- 
Robertson- Walker Universe is partially described by 



d^ + k 8nG 



{p + 3p) , 



(3) 
(4) 



where the dot represents a cosmic time derivative, a is the 
scale factor, p the density, p the pressure and fc is a constant 
(since we are mainly interested in the case of a fiat Uni- 
verse, we take fc = this point onwards). To fully specify 
the dynamics, one must supply an extra relation between 
pressure and density, in other words, an equation of state 
p = p{p) = w{p)p. 

Since the cosmological principle is only approximately 
valid, so is the description given above. It applies only to 
the average Universe. In (|3} and Q, p and p refer to av- 
erage global measures of pressure and density obtained by 
smoothing over large enough regions, as opposed to local 
measures. Thus, it would probably be better to write them 
as (p) and (p). Note, however, that the extra prescription 
p — wp is a relation between the local pressure and density 
and not necessarily between the smoothed versions. This is 
a subtle and important point. 

Consider the perturbative decomposition of pressure 
and density 



P = (p) + Sp + ■ 
p = {p) + Sp + ^ 



(5) 



The common procedure is to assume that (p) = w{{p)){p). 
However, a priori, there is no reason to suppose that the re- 
lation between the average pressure and the average density 
is the same as the local one (for a related discussion see Ellis 
& Buchert (2005). The Chaplygin gas is an example of this. 
While locally it behaves as p = —Ap~°', globally it does not 
(unless a = Q): 



(P) 



-A{p-") ^-A<p>- 



(6) 



except if perturbations are linear {S — Sp/{p) <C 1) in which 
case we would have 



(P) 



-A{p- 

'A{pr 



') = -A{p)-'^{l-a5) 



(7) 



This discrepancy between local and global equations of state 
depends on what is happening to the fiuid on small scales 
meaning that non-linearities are deeply involved. We fur- 
ther discuss this in Section 7. The important point is that 
this discrepancy complicates matters substantially. For in- 
stance, in the case of the Chaplygin gas, most researchers 
have naively used the local equation of state to relate (p) 
and (p) which in general is clearly not a valid assumption. 



3 THE SIMPLEST UDE MODEL 

The simplest unified dark energy model one can possibly 
conceive is that of a perfect fiuid with a constant negative 
pressure: p = —A with constant A > 0. We can also state 
this by saying that w = —Ap~^. Incidentally, this is the 
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case of the Chaplygin gas with a 
is straightforward to show that 



p + 3-{p + p) 



0. 



0. From ^ and Q it 



(8) 



This has a simple solution if p is constant. In fact, we have 
that p + p cc and ((3]) can thus be rewritten as 

tf = H^(n^*a-3 + ni*), (9) 

where H = d/a is the Hubble parameter and = l — O!^ = 
87rGp/3Ho is the equivalent A energy fraction today. It is 
also easy to show that wo = (p/p)o ~ — f^A* today. Here the 
index '0' indicates that the variables are to be evaluated at 
the present time, to- 

The most important feature of this model is that it is 
totally equivalent to ACDM to all orders (Avelino, Bega, 
de Carvalho & Martins 2003). This fact translates the three 
year WMAP constraint on Q'i, = 0.238t° o4? (Spergel et al. 
2006) into wo = -fit = -0.76215;;°^°. We note that the 
evolution of the equation of state of this UDE fluid is given 
by 



wo 



(1 + wo)a ^ — Wo 



(10) 



making it indeed behave as CDM {w ~ 0) at early times 
(a ~ 0) and as a cosmological constant much later (a — > oo). 



4 SOUND SPEED AND BACKGROUND 
DYNAMICS 

The simplest UDE model just described has a null square 
sound speed Cs = but that will cease to be the case in 
the context of more general models. Indeed, it is a simple 
matter to show, using Q and (U) that 



P 
P 



3 ^dH 

AttG dt 
1 d 



AtvG dt 



f'd 1,,' 
\a 2 



implying a sound speed of 



dp Id 
dp " 3HdH 



^-2 



(11) 
(12) 

(13) 



where q = — a/(aH^) is the usual deceleration parameter. 
Let us assume that H is a decreasing function of cosmic 
time. The sign of the sound speed squared today will 
be determined by the way q is evolving. If it is evolving 
sufficiently fast (towards negative values) > 0; otherwise 
< 0. On the other hand, the evolution of q is linked to how 
fast the transition from dark matter to dark energy occurs 
for quartessence. If it is steep enough (faster than in ACDM, 
that is), will be positive (negative otherwise). Therefore, 
the sign of is connected to the background dynamics. 



5 EXTENDING THE SIMPLEST UDE MODEL 

A straightforward generalization of the simplest UDE model 
leads to w = —Ap'^^^"^ with a constant a > 0, i.e. the 
generalized Chaplygin gas. It is a simple matter to show 
that the sound speed square of this fluid is given by 



dp 
dp 



(14) 



and is, in fact, positive for a > 0. This means that the 
transition from dark matter to dark energy is fast enough 
to make > 0. As we mentioned before, one interesting 
characteristic of the Chaplygin gas is that it has a non-zero 
minimal density. No matter how much you expand it, p will 
never drop below a certain value: pm = When 
p reaches this value, however, the pressure will be exactly 
—pm, that is, w = —1; it behaves as a cosmological constant. 



6 ONE-TO-ONE MAP BETWEEN UDE AND 
QUINTESSENCE 

Let us now consider the simplest quintessence model with 
PQ ~ wqpq with constant wq plus pressureless CDM and 
make a one-to-one map into a unifled dark energy model 
with density p = PQ + Pcdm and pressure p = PQ+Pcdm = PQ 
(see also Kamenshchik, Moschella & Pasquier 2001; Berto- 
lami. Sen, Sen & Silva 2004). If the density perturbations 
are small then the sound speed of this fluid is uniform and 
is given by 

^2 ^ dp _ dpQ dpQ 
dp dpQ dp 

1 

= WQ- 



(15) 



"Q 



. , (K/(l + mQ))a3'"c 

where k = p'idml Pq ~ 0.3 — 0.4. The equation of state of 
this fluid has the following simple form, 

\ i/(i+»g) 

P=—+ Pcdn, I „ ' I ■ (16) 



P 



K^QPq 



If we the dominant energy condition {w = p/p > —1) to be 
valid at all times then we must have wq > —1. In this case 
< at all times. 



7 NON-LINEAR EFFECTS 

General relativity is a non-linear theory of gravity which 
means that, unlike in the electromagnetic case, the super- 
position principle is not valid. However, in many situations 
we may neglect non-linear terms and effectively linearize the 
field equations. In cosmology, this is often possible for large 
enough scales. However, this is in general not the case in the 
context of perfect fluid UDE models. 



7.1 Case with cf > 

In this section we illustrate in the simplest possible manner 
how non-linearities can affect the evolution of the Universe 
even on large cosmological scales, focusing on the particu- 
lar example of the Chaplygin gas. Here, one should bear in 
mind than the Chaplygin gas has a minimum density pm at- 
tainable and consequently a minimum pressure pm = —pm 
(corresponding to a maximum in modulus). 

Consider a spherical region of radius R with an average 
density (p). If the energy density is uniformly distributed 
then we clearly have (p) = —A{p)~°' . However, if the density 
is not uniformly distributed then, in general, this will no 
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longer be valid. Take the case in which the region with _Ri < 
r < R has a smaller density than the region with r < Ri . In 
fact, let us assume that p(r < Ri) = N{p), where is some 
factor higher than 1, and p{Ri < r < R) = pm = A^/^^+"^ 
so that p{Ri < r < R) = pm = —Pm- The sum of the masses 
inside the two regions divided by the entire volume still has 
to be (p) by construction which implies that N is equal to 




Now, while the average density inside R is still the same £is 
in the uniform case, the average pressure (p) is not. It is 
simple enough to show that 

where the approximation is valid for large A'^ (small Ri/R). 
If N ^ 1 , the real average pressure (p) , is considerably larger 
(in modulus) than A{p)~'^, the modulus of the average pres- 
sure of the Chaplygin gas in the absence of perturbations. 
Recall that for the Chaplygin gas, the lower the density, the 
bigger the pressure (in modulus). This is why the pressure in 
the lower density region will dominate the average pressure 
inside R. 

Hence, the non-linear collapse will make (p) = pm = 
—pm early on, and will thus anticipate and slowdown the 
transition from dark matter to dark energy, thus mimicking 
the background evolution in the simplest UDE model with 
a = (or, equivalently, the ACDM scenario). The small 
scale-structure of the Chaplygin gas is interfering with the 
evolution of the Universe on large scales. 

Of course, this is an oversimplified picture: we have not 
taken into account the dynamical effects of pressure gradi- 
ents. In high density regions the pressure will be significantly 
smaller (in modulus) than average. Still, we need to take into 
account that gravitational collapse will only be effective on 
a given scale A if A > CsH~^ . However, note that once the 
perturbations become non-linear pressure will not influence 
significantly the subsequent dynamics. 

7.2 Case with cl <Q 

Let us recall what a negative cj means; remember that well 
inside the horizon, linear theory describes the evolution of a 
perturbation as a wave: 

5-c^V^<5~0, (19) 

so that 5k oc exp {i{wt — k • r)) where = c^k^ (although 
this is a linear result, we will still use it to guide us into the 
non- linear realm) . This makes the interpretation of the sign 
of the sound speed very straightforward. If > 0, 5 will 
oscillate as an acoustical wave, acting against the formation 
of voids and dense regions. On the other hand, if < 0, 
the opposite will happen: density perturbations in collapsing 
regions and voids get amplified. 

Let us now apply to the UDE model of Section 6 the 
same reasoning we did to the non-uniform Chaplygin gas 
and focus on a background of negative < (the relevant 
case) . High density regions will tend to behave as pressure- 
less matter with ~ 0. On the other hand, low density re- 
gions with < will tend to become increasingly emptier. 



There is, however, a major difference to the Chaplygin gas 
case: there is not a positive minimum density underdense re- 
gions cannot go below (except if w = —1) and, consequently, 
p can be arbitrarily close to zero. This has the interesting 
consequence that the average pressure may be close to zero 
at all times so that the Universe may never start to acceler- 
ate (despite the linear theory prediction), something which 
is clearly inconsistent with current observational evidence. 



8 CONCLUSIONS 

In this paper we have shown that there is a one-to-one 
map between any dark energy model and a perfect fluid 
quartessence model. Note, however, that in general these 
models do not share the same underlying physics. Also, al- 
though these models are equivalent in the absence of pertur- 
bations, this will in general not be the case in the presence of 
density fluctuations. In particular, if the perturbations are 
small then the evolution of cl for the quartessence model 
is fully coupled to the large-scale dynamics of the universe. 
However, we have shown that non-linear effects are of capital 
importance and should be included in quantitative treat- 
ments of perfect fiuid UDE. In fact, small scale structure 
may alter the background evolution of the Universe making 
the global behaviour of quartessence quite different from the 
one predicted by linear theory. In particular, we have found 
that if Cg > then the transition from the dark matter to 
the dark energy dominated epoch may be anticipated with 
respect to linear expectations, leading to a background evo- 
lution similar to that of standard ACDM models. On the 
other hand, if < the transition from the decelerating 
to the accelerating phase may never happen. There is one 
obvious consequence of these results: if future observations 
turn out to be incompatible with the background evolution 
predicted in the context standard ACDM models then that 
will also constitute a serious challenge to perfect fluid UDE 
scenarios. 
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